Impact of optical phonon scattering on magnet otransport in 
double-barrier heterostructures 

Dae Kwan Kim* and Patrick Roblin^ 

Department of Electrical Engineering, The Ohio State University, Columbus, Ohio 43210 

Kwang-Sup Soli* 

Department of Physics, Seoul National University, 
Seoul 151-742, Korea 

Chul Koo Kim§ 

Department of Physics and Institute for Mathematical Sciences, Yonsei University, Seoul 

120-749, Korea 

Abstract 

The transport in double-barrier heterostructures of electrons interacting with 
longitudinal optical phonons in the presence of parallel electric and mag- 
netic fields is analyzed theoretically with the aid of a 3-dimensional quantum 
transport simulator. Inter Landau state transitions induced by LO-phonon 
scattering-assisted resonant-tunneling is shown to be an important process 
with a probability comparable to that of intra Landau state scattering. Anal- 
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ysis of the current-voltage characteristics reveals also that the current peak is 
a periodic function of the inverse of the magnetic field, with a period depen- 
dent on the quasi-resonant energy level. 
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I. INTRODUCTION 



The magnetotransport phenomen, 



a0 in multiple-barrier heterostruct 



ures in which the 



magnetic is applied parallel to the electric field, was first reported by Mendez, Esaki and 
Wang§. The experimental data acquired in the presence of a magnetic field exhibited ab- 
normal features such as a shift of the onset voltagei, an increase of the peak-to-valley 
current ratio, as well as the presence of additional shoulders in the valley region of the I-V 
characteristics^. However, these observations have not been fully accounted for on a theo- 
retical basis. This is partly due to the fact that a 3-dimensional treatment of the electron 
motion is required for rigorously analyzing the impact the quantization in the transverse 
plane has upon the longitudinal electron motion. Indeed relatively few works analyzing het- 
erostructure devices with a magnetic field parallel to the electric field have been reported 
until nowBi. 

When the magnetic field and electric field are perpendicular to the interface plane of 
the composing heterostructure materials, the density of state (DOS) of the electrons in the 
plane perpendicular to the superlattice direction assumes discrete values corresponding to 



and this strongly impacts the electron motion. The resulting transmission probability and 
current-voltage characteristics of the heterostructure studied become then quite different 
from the case with zero magnetic field. 

The additional shoulders observed in the valley region of the I-V characteristics in 
the presence of a mag netic fieldi'li have been attributed to the longitudinal optical (LO) 
phonons which is one of the dominant scattering processes of electrons in the III-V semi- 
conductor heterostructure. We present here a model of resonant-tunneling in the presence 
of of a magnetic field, which includes the electron-LO phonon interaction!, so as to account 
for the transition of the incident electron to different Landau states. All possible transi- 
tions, including intra-Landau state scattering will be considered in our model and numerical 
analysis. 




When the magnetic field is varied, the DOS of the electrons changes 
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The effect of the magnetic field is to quantize the energy levels of the electron into Lan- 
dau levels. This is due to the fact that the projection of the helicoidal trajectory of the 
electron on the perpendicular plane corresponds to an harmonic oscillator. To account for 
the effect of this quantization on the longitudinal motion of electrons in the superlattice 
direction, we will develop a 3-dimensional transport model. We will then study the abnor- 
mal magnetotransport features with the aid of a 3-dimensional quantum simulator which 
numerically accounts for the contribution of the electron states in the perpendicular plane0. 

The heterostructure device considered consists of an undoped Alo.3Gao.7As / GaAs / 
Alo.3Gao.7As, double barrier structure sandwiched by heavily doped n + -GaAs left and right 
contact layers. Our analysis only considers tunneling from electrons in the T valley. The 
T conduction-band edge in the absence of an applied voltage for the device studied is 
shown in the inset of Fig. |l| Spin splitting in the magnetic field is ignored for simplic- 
ity. Finally in this paper we are interested in the region of operation for which we have 
hw c (cyclotron frequency) ^> kT, and our numerical calculations will be consequently per- 
formed at the low temperature of T = 4.2K and for magnetic fields up to 22 Tesla. 

In section II, the 3D magnetotransport model used for electrons interacting with longi- 
tudinal optical phonons via polor scattering is introduced. The electron state is expanded 
in the generalized Wannier-Landau basis. Expression for calculating the device transmis- 
sion coefficient and current-voltage characteristics in the presence of phonon scattering are 
derived. In section III, the numerical results obtained for the main abnormal magnetotrans- 
port features using the above theoretical model are presented and analyzed. Finally in the 
last section, a brief summary of our results is given. 



II. THE TRANSPORT THEORY 



4 



A. Hamiltonian 



The Hamiltonian of an electron system interacting with longitudinal optical phonons 
under a magnetic field in the double barrier structure is written as 

H = H e + H e ^ph, (1) 

where the free and interaction terms are given by 

# e = £e(A)4c A , 

A 

H e . ph = * £ a q {a q e^ r - a{e^}, (2) 

V * ' q 

where c\, c\ (a q , a q ) are respectively the destruction and creation operators for electrons 
(phonons) and a q is the interaction weight for the phonon wave vector q. Here A stands for 
the longitudinal crystal wavevector k x and Landau level L and k y in the perpendicular plane, 
which results from the choice of Landau gauge! for the vector potential A(r) and q is the 
mode of longitudinal optical phonon. The LO phonons will be represented by an Einstein 
model with the constant phonon frequency w. The Schrodinger equation to be solved is 

1 1 dt 1 ' w 

We expand the electron wave function in terms of generalized Wannier functions in the 
superlattice direction, x-direction0'0 and Landau states in the (y, z) plane, 

I*) = E/ dkyf(k y ,L,n)\k y ,L,n) (4) 

n,L 

To solve the Schrodinger equation, we calculate the matrix elements of the Hamiltonian 
in the generalized Wannier-Landau basis. For the free part of Hamiltonian, we find 

(k y ,L,n\H e \k' y ,L',n') = H* nn ,(L)6 L ,, L 5(k y - k y ), (5) 

where 

H e nn ,(L) = H nn , + (L + l/2K(n)<W, (6) 



with 



H nn > i w n+l,n' 



1 

2d 2 \J 'm* (n)m* (n + 1) 

h 2 

2d 2 y l m* {n)m* {n — 1) 

h 2 

— — -y-r + E con (n) - eV app {n) 



S n -l,n' 

$n,n> (7) 



a 2 m*(n) 

and w c (n) = TieB / m* (n) , the cyclotron frequency at the lattice site n. Here, E con (n), V app (n) 
are the conduction band edge and applied voltage at the lattice site n, respectively. The 
band model selected consist of a tight-binding band structure in the superlattice direction 
and a parabolic band structure in the perpendicular plane. 

Note that the transversal mass in the Hamiltonian matrix varies in the longitudinal 
direction. It is convenient to rewrite the Hamiltonian matrix element as 

H e nn ,(L) = H nn ,{L) + (L + 1/2K(0)<W, (8) 

in which the position (n) dependent transversal component is now absorbed in 

H nn ,{L) = H nn , + (L + l/2)hw c (0) - l) 5 nn ,. (9) 

\m*(n) J 

H nn i which depends on the perpendicular states L, incorporates the 3-dimensional effects in 
the longitudinal motion of the electrons. 

After some algebra, the matrix element of the interaction Hamiltonian which determines 
the inter Landau state transitions leads to 



1 E 1 

q x >o 

^{a^ y (-q ± )e^ t+ M + a + F K L f L y (+ q± )e^ wt+ ^}, (10) 



(k' y ,L',n'\H e - ph \k y ,L,n) = — = J2 ~ sm(q x na)5 n/n 

q± 

where w is the phonon frequency, q± is the transversal part of q, and q is the phase of the 
LO phonon of mode q (for a more detailed formalism, see ref. [9]). The coefficients a± is 
given by 
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a_ = (2eF/e Q )y/N ph + l and a+ = (2eF / e ) ^/ N ph , (11) 

which corresponds to the emission and absorption of phonons, respectively. Here, eo is the 
vacuum dielectric constant, and N p h is the equilibrium number of phonons. F is given as 
F = (hw LO el/2) ■ {l/e opt - l/e stat ) with phonon energy hw LO = 35.3 meV and e sta t/e opt = 
1.1664 for GaAs. At liquid Helium temperature, the number of optical phonons N p h is very 
small, so only the emission of phonons is possible. In Eqn.(lO), we introduced the coefficients 

F k jt(±q±) = (k' y ,L'\e^\k y ,L) 

= 5(k' y -k y ± gj/ ) e [±^(fcL-*»±9»)i 2 /2+i(iV-iV')(flW4)] j p i ., L ( 9 2 /2/2)) (12) 

with 

where N = max(L, V) and N' = min(L, V) are the maximum and the minimum of L and 
L', respectively, 6 = sign(N — N') arctan(q y /q z ), l 2 = h/eB, and L^~ N (x) is an associated 
Laguerre polynomial!!. 

B. The envelope equation 

Using the obtained matrix elements, Eqn.(5) and (10), in the above subsection, 
Schrodinger equation can be reduced to the following 3-dimensional envelope equation for 
f(k y ,L,n), 

ih^f(k y , L, n) = (L + l/2)hw d0 f(k y , L,n) + Y^ H nm f(k y , L, m) 



i Qx>0,q ± i 

— y - 



a.^{q x na)f{k^L',n)Fl^{-q ± )e^ t+ ^ 

(14) 



+ a + sm(q x na)f(k' y , L', n)F k L y fi (+ qi _)e^ t+ ^ 
where w c q = eB/m*(0), is the cyclotron frequency at the left contact. 
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Let us label the quantum state of an incident electron injected at the left contact with 
the energy E = E 0x + (L + l/2)hw c0 using the quantum numbers (E 0x , L , k y0 ). An electron 
in the incident ballistic state (Eq x , Lq, k y o) which is scattered by a phonon will transit to 
other accessible states (E x , L, k y ) with total energy E = E ± Tiwlo following the emission 
(-) or absorption (+) of a phonon. Therefore the solution of the envelope Eqn.(|TJ]) is chosen 
to be of the following form 

-iE t 



f(k y ,L } n) = e-^ ot f (n)5(k y - k y0 )5 



LL 

+ ^=e-W)< £ -f + (n,q±,q x )F k L l k ;°(q±)e-^ 



+ ^=e-^-^ Yl -Un,q ± ,q x )F k L l k : (-q^, (15) 



" " ^ -J [<>.q..q,)r' 

where w is the phonon frequency w = wlo, fo is the coherent solution for the envelope 
equation and where f± are the scattered waves due to absorption or emission of a phonon, 
respectively. After substituting this expression into the envelope equation and taking an 
ensemble average0, we arrive at 

E 0x f (n) = Y J H nm f (m) + J2 { h se( u i m ) + h se( u i m)}/oH, (16) 

m q x >0,m 

f±(n, E lx , q x ) = Y G nm (E lx ) sm(q x ma) f (m) , (17) 

m 

where H SE is the self energy due to absorption and emission of phonons. The self energy is 
verified to be given by 

a 2 q x >0,q ± ^ 

Hg E (n,m) = -± Y —sm(q x na)sin(q x ma)G nm (E lx )[E L(hU (q 2 J 2 /2)} 2 , (18) 
si jj q 

and G nm (Ei x ) is the longitudinal Green's function in the generalized Wannier basis evaluated 
using the relation 

z2(Ei x — H) n iGl m = 5 nm . (19) 
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C. Transmission probability and current density 



For an incident electron state (E 0x ,L ) the current entering on the left hand side and 
leaving the quantum region on the right hand side are respectively given by 

j I0 (E 0x ,L ) = ev L (E 0x ) 

j T (E 0x ,L ) = e\f (n R )\ 2 v R (E 0x ) + j±, ph (E 0x , L ). (20) 
Here the phonon-assisted current component J± iP h is given by 



2 Qx>0,q s 

j±* h (E 0x ,L ) = e-£ £ ^\f ± (n R ,q ± ,q x )\ 2 [F Lo:L/ (qll 2 /2)] 2 . (21) 
sZ j^i q 



In the above equations, Vl and v R are the electron velocity at the left and right contact, and 
tir is the lattice-site index of the right contact. The transmission probability from the left 
to right contact is then defined as 

rr ,-ns t \ [jTo(Eox, Lq) + j +!ph (E 0x , Lq) + j-, p h(E x, L )} 

1l^r{Eo x ,L ) = — — — (22J 

Before calculating the current, the electron chemical potential n(B) at the left contact 
corresponding to the doping density n in presence of the magnetic field B, should be de- 
termined. As is shown in Fig. [I], the chemical potential is found to exhibit an oscillatory 
behavior in 1/5 as the magnetic field varies. The chemical potential fi(B) in the left contact 
is obtained from the carrier density through the following integration relation, 

n = N D = (V2^/n%)D(B) £ f M^^B, L )) ^ ^ 

l q =o Jo V Eox 

where D = eB/h is a degeneracy factor in the plane perpendicular to the superlattice 
direction and where the Fermi-Dirac function Jfd is given by 

fFD(E 0x ,{i e f f (B,L )) = -— — y—p—z (24) 

exp[(E 0x - He ff {B, L Q ))/kT\ + 1 

Here we have found it convenient to introduce the effective chemical potential fi e ff for 
the Landau sub-band Lq which is related to the actual electron chemical potential n{B) 
according to 



ti ef f(B, L ) = fi(B) - (L + l/2)hw c0 (B). (25) 

fi e ff can be interpreted as the effective chemical potential of the electrons in the Landau 
state L in the presence of a magnetic field. Thus, the total current density is 

00 roc 

I tot = (eD(B)/nh) ]T / T L ^ R (E 0x , L )f FD (E 0x , fx eff (B, L ))dE 0x -[£<-> R], (26) 

where the second term in the right-hand side is the backward contribution to the total 
current. At very low temperatures such as 4.2 K, the current for a occupied sub-band Lq 
reduces to 

hot,L « {eD/irh) / dE 0x T L ^ R (E 0x , L ) -[L++R). (27) 
Jo 

Therefore, the area from to fi e ff of the transmission probability determines the contribution 
of each occupied sub-band to the total current. 

III. NUMERICAL RESULTS 

As mentioned above the heterostructure device considered consists of an undoped 
Alo.3Gao.7As/GaAs/Alo.3Gao.7As, double barrier structure (see inset of Fig. [I]) sandwiched 
by heavily doped n + -GaAs left and right contact layers. The barriers and well are both 50A 
wide. The donor density in the contact regions is 1 x 10 18 /cm 3 . The donors are assumed to 
be completely ionized due to the heavy doping density. The value of the electron effective 
mass used for GaAs and Alo.3Gao.7As is m* = 0.067mo and 0.092m , respectively. 

The numerical calculation starts with the calculation of the matrix element Fl'l(q±1 2 /2) 
which is required to evaluate the self energy H^ E (n,m) associated with the inter Landau 
state transitions for all possible values of (Eq x , Lq). This self energy is then used to solve the 
envelope equation, and obtain the transmission coefficients. The current density is finally 
obtained by integrating the transmission coefficient obtained over the incident energy. We 
confined our numerical analysis to a single-sequential phonon scattering event since this is 
a quite reasonable approximation for the small heterostructure device under investigation. 
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In order to calculate the currents, it is required to calculate the chemical potential fi 
as a function of the applied magnetic field B. Using Eqn. (p3|) the chemical potential fi(B) 
calculated for several temperatures is plotted in Fig. |1] versus B. The chemical potential 
fi(B) is found to be a periodic function in 1/B. This figure clearly indicates that in the 
region of strong magnetic fields (B > 10 Telsa) using the approximation of a constant Fermi 
energy Ef would introduce a rather large error in the chemical potential. This would in 
turn adversely affect the accuracy of the device current calculated as is verified in Fig. £|. 
However at the temperature of 100-fT, the oscillating behavior in /i disappears completely 
due to thermal population of all the Landau states near the Fermi energy. 

In Fig. 0, the transmission probability in a flat band (no applied voltage) is plotted for 
two magnetic fields of B = 10 and 20 Telsa at T = 4.2K. At this low temperature, the 
emission of phonons dominates since the number of phonons in the system is negligible. The 
peaks beyond the main peak centered at 94 meV, are phonon-assisted resonant tunneling 
peaks which satisfy the following energy conservation relation 

E 0x — hw = E x + AL ■ hw c0 , (28) 

where (E 0x ,L ) is the incident electron state and (E X ,L) is the scattered state after the 
emission of one phonon and with AL = (L — Lq). The energy of the phonon is chosen 
to be hw = 35.3 meV for GaAs contact layers for simplicity. A single-sequential phonon- 
emission has been assumed. This was reported to be a valid approximation in a small size 
system like the double barrier structure considered!!! The peak indicated by AL = in 
Fig. arises when an incident electron moves in the structure with the emission of a single 
phonon and the Landau level remains conserved while the peaks AL = 1 and 2 indicate 
that phonon-assisted resonant-tunneling induced a transition to higher Landau levels. As 
the magnetic field is increased, the spacing hw c (B) between AL neighboring peaks increases 
in proportion to the applied external magnetic field where as the spacing between the main 
peak and AL = peak remains constant as it is set by the phonon energy hw. For an 
incident state with Landau levels Lq > 0, transitions of the type AL = —n arise (with n 
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an integer verifying n < Lq), and additional peaks appears between the main and AL = 
peaks. 

The matrix element FiiL^q^l 2 /2) between Landau states in Eqn.(13) determines the 
strength of the inter Landau state transitions in the transmission coefficient. For example, 
at B = 10 Telsa its argument q]l 2 /2 assume values in the range (0.025,252) depending on 
the phonon mode. After the matrix Fii^q 2 ) 2 /2) is summed over all possible phonon modes 
q±, the inter Landau state transition probability is obtained. The analysis of the various 
scattering-assisted resonant-tunneling processes induced by LO phonon emission, indicates 
that the transmission probabilities for scattering events involving the transition between two 
different Landau states (AL > 0) is of the same order as the transmission probability for 
intra Landau state (AL = 0) scattering events. Fig. [| clearly shows this property. 

Fig. ^ f|, and || shows the electron current density calculated under various conditions. 
The current density is obtained from the integration over the incident electron energy of the 
transmission coefficient T(E 0x , L ) multiplied by the Fermi-Dirac occupation function which 
is almost a step function at the low temperature considered. A broadening of the Landau 
levels up to 2 meV was considered, but there were no noticeable difference in our results for 
the calculated currents for B > 5 Telsa since we have hw c o(B = 5) = 8.6 meV ^> 2 meV. 
We observe in Fig. || and [|, the presence of small additional shoulders of almost equal width 
in the valley region of the current-voltage curve. These shoulders originate from that fact 
that the phonon-assisted resonant-tunneling peaks in the transmission coefficient are spaced 
by hw c above the main quasi- resonant energy peak E r . 

Before discussing the impact of phonon scattering, let us first briefly describe the main 
effects of the magnetic field upon the current of the resonant tunneling diode (RTD) in 
the absence of phonon scattering. As can be seen in Fig. |3| and [| the magnetic field is 
found to induce a shift of the onset voltage of the RTD I-V characteristic by the amount 
AVb,on/2 = [/^(O) — /i(L?)] + hw c o(B)/2, and to enhance the value of peak current of the 
RTD. Fig. [| reveals another interesting feature regarding the dependence of the RTD current- 
voltage characteristics upon the magnetic field. In this Figure the diode current is measured 
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at Vd = 0.146 Volt which is in the current peak region of the RTD I-V curve. Note 
that phonon scattering was not included in this calculation as coherent (ballistic) tunneling 
gives the dominant contribution to the peak current. The RTD current measured at that 
fixed diode voltage is seen to exhibits an oscillatory behavior versus 1/B which increases 
in amplitude as the magnetic field is increased. The current maxima observed in Fig. |] 
for the current versus B are found to occur for values B n of the magnetic field satisfying 
the following relation fi(B n ) — {n + l/2)hw c o(B n ) = £7 r (Vo) where the quasi-resonant level 
E r (Vn) — E r (0) — Vd/2 is varying with the diode voltage. This leads to the interesting 
result that its periodicity versus 1/B which is approximately he / m* (0)[/a(0) — ^(Vb)] -1 , 
depends on the quasi-resonant energy level E r (Vo); for a smaller value of E r , one obtains 
an increased oscillatory behavior in the current-magnetic field curve. 

As can be seen Fig. ^ and [5], the effect of phonon scattering on the RTD current density 
is to shift the current- voltage characteristic downward compared to the one obtained in the 
absence of scattering. Another effect is to increase the current and to introduce additional 
shoulders in the valley region of the current- volt age characteristic. The downward shift of 
onset voltage results from the negative value of the real part of the self energy induced by 
the scattering of the electron by the phonons. For magnetic fields of B = 10 and 20 Telsa, 
the downward shifts obtained in the current- volt age characteristic are found to be 3.6 and 
4.0 mV, respectively. The small shoulders at B = 20 Telsa in the current valley region of 
Fig. |5| become larger than for the case of B = 10 Telsa. Thus in stronger magnetic fields, 
these additional shoulders in the valley region become higher and their intervals becomes 
wider, which agrees with the reported experiments!!! 

To account for the origin of the shoulder at Vd = 0.187 Volt, let us consider the inset 
of Fig. |5|. At B = 20 Telsa, only the two lowest Landau levels L = and 1 contribute 
to the diode current. The effective chemical potentials defined in Eqn.(25) for the L = 
and 1 states are respectively /i e // = 37.3 and 2.8 meV. They are indicated by two arrows 
in the inset of Fig. [5[ At Vd = 0.187 Volt, both the peaks AL = and —1 of transmission 
coefficients for the Lq = and 1 incident states are respectively included into the integration 
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range of the current Eqn.(|27|), explaining the formation of the shoulders. 

Currently, it is difficult to compare our theoretical results in a quantitative way to the 
experiments reported in Ref. 2 and 3. The reasons are as follows. First, the authors of Ref. 

2 concluded that additional experiments are needed to confirm the weak minima of their 
preliminary data in the valley current region for strong magnetic fields. Second, our test 
device doesn't have a spacer layer in the emitter region, while the device structure of Ref. 3 
have a spacer layer where a discrete triangular well state is formed. The shoulders therein 
corresponding to 5L ^ beyond the main peak was reported to be as high as almost 2/5 
of the main peak for extremely strong magnetic field, which were clarified to originate from 
the assistance of LO phonons to the electron current. Our results together with Ref. 2 and 

3 lead to the conclusion that the spacer layer tends to enhance the strength of the shoulders 
beyond the main peak. But in order to explain such high shoulders, other contributions 
due to scatterings of electrons by interface roughness, acoustic phonon, etc., might also be 
needed. 

IV. SUMMARY 

In our paper we have investigated the transport of electrons interacting with LO phonons 
in a double-barrier heterostructure with a magnetic field applied parallel to the electric field. 
The effect of the magnetic field is to quantize the energy of the electron in the perpendicular 
plane to the magnetic field. To account for the effects of a magnetic field applied in parallel 
with the superlattice direction, a 3-dimensional transport model was developed which ac- 
counts for the variation of the transverse mass and the associated variation of the cyclotron 
frequency across the heterostructure. 

In order to determine the inter Landau state transitions in the transmission probability, 
all the matrices -Fzy^g^Z 2 /2) for allowed phonon modes have been calculated numerically. 
It was established that the inter Landau state transitions due to LO phonon scattering are 
as probable as intra Landau state scattering. 
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For a fixed diode voltage in the peak region of the RTD current- volt age curve, the RTD 
current versus the magnetic field was found to be a periodic function of 1/B. The current 
maxima have an 1/B period which is approximately given by he/m*(0)/(fi(0) — EriVo))- 
The dependence of this 1/B period on the quasi-resonant energy level £V(Vd) of our test 
double barrier structure indicates that magnetotransport can be used for the spectroscopic 
analysis of RTDs for diode voltages in the peak current region. 

A theory and associated numerical analysis using scattering matrix elements was pre- 
sented to describe the impact of electron-LO phonon interaction in the I-V curve in the 
presence of both parallel electric and magnetic fields. The results obtained were found to be 
in qualitative agreement with the experimental results reported by Ref. 2 and 3. However 
for further understanding and improved quantitative fit of the weak shoulders in the valley 
region of the double barrier IV, more theoretical work on magnetotransport including scat- 
tering processes such as interface roughness and acoustic phonon scattering, etc., is required 
in addition to LO phonon scattering studied in this paper. 
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FIGURES 

FIG. 1. The chemical potentials at n + = 1 x 10 18 /cm 3 as a function of the magnetic field for 
various temperatures. The solid, dotted, dashed, dot-dashed, double dot-dashed lines correspond 
to T=4.2, 10, 20, 40, and 100 K, respectively. The bold dotted line is the Fermi energy Ef at 
T = 0, B = 0. The inset figure is the conduction band profile of the resonant-tunneling structure. 

FIG. 2. Transmission probability for the incident Landau state Lq = with no applied voltage 
at T = 4.2K. The dotted and solid lines correspond to B =10 and 20 Telsa, respectively. 

FIG. 3. RTD current densities versus applied voltage for B = 10 Telsa at T = 4.2K. The dot- 
ted, solid lines correspond to the RTD currents calculated without and with LO phonon scattering 
respectively, and the dashed line is for the RTD current in the absence of both an applied magnetic 
field and phonon scattering. In inset a plot showing the I — V characteristic in the absence of both 
an applied magnetic field and scattering for a wider diode voltage range is also given as a reference. 

FIG. 4. Current density (solid line) plotted versus the magnetic field for an applied voltage of 
0.146 Volt. The dotted line corresponds to the RTD current density calculated in the approximation 
of a constant Fermi energy (//(-B) = m(0)). The dashed line in the center corresponds to the current 
value in the absence of a magnetic field (B = 0). 

FIG. 5. RTD current densities versus applied voltage for B = 20 Telsa at T = 4.2K. The dot- 
ted, solid lines correspond to the RTD currents calculated without and with LO phonon scattering, 
respectively. The dashed line is for the RTD current in magnetic free case without scatterings. 
The inset shows the transmission probability for the two lowest Landau states Lq = and 1 at the 
voltage V D = 0.187 V. 
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